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Abstract. The exact dynamics of a disordered spin star system, describing a cen- 
tral spin coupled to N distinguishable and non interacting spins |, is reported. 
Exploiting their interaction with the central single spin system, we present possible 
conditional schemes for the generation of W-like states, as well as of well-defined 
angular momentum states, of the N uncoupled spins. We provide in addition a 
way to estimate the coupling intensity between each of the N spins and the central 
one. Finally the feasibility of our procedure is briefly discussed. 

1 Introduction 

Interacting spin models play a central role in many physical contexts providing a paradigm to 
describe a wide range of different systems. In condensed matter physics, for example, they can 
be explored to analyze many properties of magnetic compounds. It is indeed well known that 
a suitable general model of a magnet consists of N spins coupled by exchange interaction with 
arbitrary range and strength. 

The interest toward spin systems, and more in particular toward spin dynamics in semicon- 
ductor structures, has remarkably increased in the last few years also in connection with the 
new emerging areas of quantum computation and information pQ-[3]. In these contexts, spin 
models like Heisenberg spin chains or spin star systems, describing for example a single electron 
spin in a semiconductor quantum dot interacting with surrounding nuclear spins via hyperfine 
coupling mechanisms, have been extensively studied [4]-[14j. 

Generally speaking spin models have proved to be promising candidates for the generation 
and the control of assigned quantum correlations, as witnessed by the numerous papers recently 
appeared in literature [15] - [18] . 

In this paper we concentrate on the possibility of manipulating at demand the state of a 
sample of N uncoupled spins exploiting their common interaction with another single spin called 
central system. The analysis we have developed, on the one hand provides possible procedures 
to guide the system toward pure states characterized by fixed correlation conditions, on the 
other hand suggests a way to to estimate the coupling strenght between each of the N spins 
and the central one. 



2 Disordered spin star system 

Having in mind as objective the possibility of generating fixed entanglement conditions in a 
system of N uncoupled spins i> we exploit the interaction of each element of the system with 
another spin i hereafter called central one. The system we are talking about is illustrated for 
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Fig. 1. An illustration of a disordered spin star system 



convenience in Figure 1 and it can be described adopting the following hamiltonian model: 

H = H + H! (1) 

with 

JV 

H = u)J2 a i+"oVz, (2) 
i=i 

JV 

H I= J2^(aial+a A ai), (3) 

3=1 

The Pauli operators labelled by the index A refer to the central spin, while the others char- 
acterized by the index j (j = 1..N) refer to the N spins. The coupling strength between the 
spin j and the central one is measured by the constant aj and, generally speaking, ctj may be 
different from for j ^ i. In realistic physical situations the coupling constants can change 
for example proportionally to the distance from the central system . We refer to the system 
described by eq. ([T]) as disordered spin star system. 

Studying the dynamical properties of this model is for example of interest in contexts like 
quantum dot coupled by hyperfine interaction with nuclear spins, or electronic spins bound to 
phosphorus atoms in a matrix of silica or germanio in presence of defects |20) . 

The symmetry properties of the hamiltonian model can be successfully exploited in order to 
analyze the dynamics of the system. It is easy to convince oneself that the component along the 

(T A 1 JV " <J A 

z axes of the total angular momentum operator S z = -j- + ^ X^=i = -|- + J 2 is a constant 
of motion. Then, starting from an eigenstate of S z the system evolves in the correspondent 
invariant Hilbert subspace. Let's suppose in particular to prepare the system in the following 
state 

\Tp{i 1 ,... ip} (o)) = \n)\i..-u 1 ..-u P -..i), (4) 

where the central spin, as well as p of the N spins, namely the ^-th, i 2 -th,... i p th, are in their 
respective up state |f) defined as <j\ |t) = |T) with i = A, 1, .., AT, whereas the others are in their 
down state \[) with a\\[) = -||). 

We in addition denote the state |Ta)|I, li • ■ ■ > !)i where all the uncoupled spins are down, 

by|Vo(0)>. 
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Taking into account the previous considerations we may claim that, at time instant t the 
state of the system prepared in the state ((4]) can be written as 

^{i!,.. .*„}(*)) = y. j: . ,,. .... ,. ajuj2,...,j P (t)\U)\l ■ ■ • Tii ■ • • hp •••!> + 

+ Sj 1 <j 2 <...<j p+1 bj 1 ,j 2 ,...,j p+1 {t)\lA)\l ••• Tii ■ ■ • I ■ ■ • Tip+i •■■I), (5) 

In equation ([5]) each index ji (i — 1, ..,p + 1) runs from 1 to N. Thus the first term of the 
right hand side is a superposition of all the states in which the central spin, as well as p among 
the N uncoupled spins, are in their up state. The second term of equation |5j) is similarly a 
linear combination of all the states in which the central system is in its down state whereas 
p+ 1 of the N spins are in their up state. We notice that the structure of \ipo(t)) can be deduced 
from that of |V'{ii....i p }( i )) simply substituting to the first sum of eq. (fSJ) the term a(t)\^A) 
|| ... 1 ... I). Inserting eq. ([5]) in the time-dependent Schrodinger equation leads to the following 
system of coupled equations for the probability amplitudes aj 1 ,j a ,...,j p (t) and bj 1 j 2i ....j p+1 (t) 

N 

ih juh,---,jp( t ) = Aa h,h,---,jp( t )+ a r b o({h,...,j p }u{r})(t) (6) 

N 

X] a r a Mii,...jp+i)W ( ? ) 

r=l(re{ii,..i„}) 



^ii.i2,...,ip+i W — A b ji,h,---,jp+i W + 



where A = uj — uj . In eq.© we have introduced the operator O which adds the index r to the 
set of indices {ji, . . . , j p } arranging them in increasing order. We point out that this operator 
is well defined if r does not belong to the set {ji, . . . ,j p }, otherwise the probability amplitude 
bo({ji,...,j p }u{r})(t) would have p indices instead of (p+ 1) becoming senseless. The operator 5 r 
appearing in turn in eq. (|7j) acts on the family of p+ 1 indices, recovering a set of p indices from 
{jii J2, • • ■ j jp+i} by eliminating the index r. We have to mention that the above operator is 
well defined if r belongs to the set {ji, J2, ■ • • , jp+i} m order to assure the correct definition of 
a probability amplitudes of the type a>j 1 ,j 2 ,...,j p (t)- 

This system of differential equations can be easily decoupled when the system is prepared 
in the state l^o (()))■ In this case eqs. (j6|) and ([7|) become 



A' 



ia(t) = Aa{t) + ^a J b ] {t) 



(8) 



ibj(t) = -Abj(t) +a 3 a(t), 
and it is easy to show that a(t) fulfills the following Cauchy problem 

a(0) = 1 



&(t) = -[A 2 +J2< 



o(t), 



(9) 

(10) 
(11) 



whose solution is 



a(t) = cos 




A 



: sin 



A 2 




(12) 



Inserting eq. (fl2")) in eq. © leads to a first order non homogeneous linear differential equation 
for bj(t) which can be easily solved from the initial conditions bj(0) — getting 



b 3 (t) 



a , 



N 



(13) 
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When the system is prepared in the state (|4|), eqs. ((6]) and (JT]) may be still managed in such a way 
to decouple the probability amplitudes of the type %'i,j 2 ,...,j» from those of type bj 1 j 3 ,...,j p+1 - 
We do not present here such a procedure since in what follows we concentrate on the rich 
dynamical properties of the system evolving in accordance with \ipo(t)). 



3 Generation of W-like states 

The results obtained in the previous Section suggest that measuring the central spin we have 
the possibility of guiding the system of interest, namely the N uncoupled spins, toward a linear 
coherent superposition of states characterized by the fact that only one spin is in the state |f) 
whereas the others are in the state |4). 

Starting from eq.([5|) we may indeed claim that a measure of the observable gives, with 
probabilities Ylj=i l^i W| 2 , the eigenvalue —1. In this case the iV spins are left in the normalized 
state 

\W^ n ) = 1 MT, |, . . . , 4) + . . . + aj\i, . . . , T, ■ • • , 4) + • • • + ajv|4, 1, •••,!)) • 

(14) 

The state given by eq. (fl~4f looks like the well known W-state [21] 

\w) = J= (|T, 4,- ..,4> + -.. + 14,-. 4) + ... + 14,4,- ••,!))• (15) 

v jV 



the only difference between the two states (|14|) and (fl5|) being the weight of each component 
in the superposition. In the W-state all the states of the superposition appear indeed with 
the same probability. On the other hand the state (|14p we have obtained with our procedure, 
reduces to the W-state when the central spin does not distinguish the N spins around it, that 
is when ctj = a Vj = 1, . . . , N. For these reasons we call the state \W gen ) a W-like state. 

It is important to underline that the procedure we have discussed is a conditional one. In 
other words we may claim to generate the state IT^ 96 ™) only if the measurement of gives the 
eigenvalue —1. Starting from eq. (|13|) we can write the probability of success of our procedure 
as follows 



N 2 
' ' X^N 2 



N 



3 



A 2 t\. (16) 



Thus, appropriately choosing the time instant t at which the measure of is performed, 
we may generate the desired state with the highest probability that coincides with one in 
correspondence to A = 0. If indeed we measure the observable at time instant 

7r(2n + l) , x 

tn= r -— ' , 71 = 0,1,2,... (17) 

the probability of success becomes P = Y", N 1 ° 3 and thus P = 1 if A = 0. 

It is on the other hand of relevance to analyze the behaviour of such a probability with 
respect to imprecisions in setting the time instant at which the measurement of the observable 

is performed. Let's first of all observe that the probability of success P, as given by equation 
(fT6|) . is a periodic function of t, the period being T — 7r(52 - =1 + A 2 )~^ . This circumstance 
suggests to choice the time instant at which to perform the measurement optimizing, as far as 
possible, our proposal on the experimental side too. To this end let's consider for simplicity the 
case A — and indicate by i* — t n + S the time instant at which the measurement if performed. 
If 6 is small enough the probability of success does not appreciably get reduced. Let's suppose 
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in particular that \S\ < T10 1 . Moreover it is also reasonable to assume that the relative error 
j- is of the order of 10 -3 thus implying that 

t„ ~ 10 3 5 < 10 2 ^ = (18) 

If our model describes, for example, hyperfine interaction of a localized electron with nuclei, 
the order of magnitude of the coupling constant ctj can be estimated as 10 _5 eV [55], [53] in 
correspondence to which t n must be less than or equal to 7rl0 _9 sec. Thus t„ ~ Ins is compatible 
with eq. (|18[) and it may be realized fixing n = 100. In figure 2 we plot the probability of success 



Ptob ability 
of success 



0.8 



0.4 




0.£ 




Fig. 2. Probability of success to generate |W /9en ) as function of x — -f- in correspondence to A — 
and n=100. 



of our scheme given by eq. (fT6|) versus x — putting A = and t = tioo + 5. As foreseeable, 
this probability of success remains of experimental interest for x up to 3 • 10~ 3 . Thus we may 
conclude that the procedure is stable enough against unavoidable uncertainties in the time 
instant at which the measurement of is done. 

Before concluding this section it is interesting to emphasize that the coupling between the N 
spins and the central system generates entanglement between any two spins around the central 
one. In order to estimate the amount of such an entanglement and analyze its time evolution 
we evaluate the relative concurrence function. It is easy to prove that, starting from the initial 
state IV'o(O)), the reduced density Pij(t) describing the system of the two spins i and j among 
the N uncoupled ones, can be written in the form 



pij{t) 









hit)] 2 bi^bjit)* 
ob z (tyb 3 (t) \bj(t)\ 2 

|a(t) | 2 



(19) 



when expressed in the standard two-spin basis {|TT)> Itl); lit); III)}- I n equation (fT9l) bj(t) and 
a(t) are given by eqs. lfT3|) and (fT2"]) respectively. Thus the concurrence function Cjj(t) turns 
out to be 



Cij(t) = 



2\aA\a, 



Z^=i 



A 2 



■ sm 



N 



A 2 t 



(20) 



As expected, the degree of entanglement get established between the two uncoupled spins i and 
j, oscillates with time t. In addition the maximum value of Cij(t) is proportional to |a,||aj| and 
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when afc = a Vfc and A — reaches the value jj. This means that, increasing the number of 
spins around the central one, leads to a weak and weak pair quantum correlation in the system 
of interest. 



4 Generation of well-defined angular momentum states of the N uncoupled 
spins 

The W-state in eq. (fTB"j) is a multipartite entangled state of the N uncoupled spins around the 
central one. On the other hand it coincides with a particular superposition of states in which 
N — 1 spins have projection down while only one has projection up. This state can be thus 
obtained applying the collective ladder operator a + = J2i a + on the state ||, |, . . . , |) that is 
|W) x Eti °+ l-lo I) •• • i-L)- I R other words the W-state is a common eigenstate of J 2 and J z , 
J = i X)j=i a i Dem g the collective angular momentum operator of the N uncoupled spins, with 

J 2 |W0 = f(f +\)\W) (21) 
J,\W) = (-% + 1)\W) (22) 

Thus \W) = | J, M) whit J = f and M = -f + 1. 

This observation suggests us the possibility to iterate our procedure in order to generate all 
the well defined angular momentum states | J = ^, M) with M = — -y . . . ^ . 

Let's indeed consider the spin star system in which the central spin interacts in the same 
way with all the others uncoupled N spins, that is ay = a Vj. Under this condition the 
hamiltonian model ([!]) is invariant by permutation of an arbitrary couple of spins among the N. 
Moreover [a A2 ,H] = [J 2 ,H] = [af + J z , H] = [jf nt ,H] = 0, 3i nt being an intermediate angular 
momentum resulting from the coupling of selected at will individual angular momentum of the 
N spins. These symmetry properties suggest to develop the dynamics of our system exploiting 
the coupled angular momentum basis { | J, M, v) } for the N spins instead of the factorized one 
previously used. The index v runs from 1 to I'm ax (J) and allows us to distinguish between 
different states of the basis characterized by the same J and M. 

In the initial condition \ipo(0)} the N spins around the central one are in the coupled angular 
momentum state |J = ^-,M — = 1)- In what follows we do not indicate anymore the 

the index v remaining it equal to one. Thanks to the symmetry properties of our system at a 
generic time instant t we can rewrite the state of the total system in the form 

\m) = Ai{t)\ U)|y,-y) + Bi(<)| U)|y, -y + 1) (23) 



where 



Ai = cos(piat), (24) 
Bi = — isinipxOit), (25) 



with 




NW N 1)=VN. (26) 



2 / V 2 

Thus as before, if we assume the central spin A measured in the state | J,^) the N uncoupled 

N _N_ 

2 ' 2 



spins are projected onto the W-state \W) = \^-,— + 1). As demonstrated in the previous 



Section the probability of success to generate the W-state coincides with sin {piott) and it 
is equal to 1 if the measurement is performed at time instants t n — 2apT^ obtained from 
equation (TT7|) putting ay = a Vj and A = 0. Suppose now to iterate the procedure preparing 
once again the central spin in the up state. The new initial condition is then 

mo)) = \u)\j-j + i) m 
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that as easily demonstrable evolves as 



with 



\m) = M{t)\ , -y + 1) + B 2 (t)\ U )\E-*L + 2) (28) 



A 2 (t) = cos(p 2 ai) (29) 
B 2 (t) = -ism(p 2 at) (30) 



where p 2 — \J y(t + -0 — ( — T + T + 2). Eq. (f2"8")) immediately implies that, measuring 
the spin A in its down state, makes the N spins to collapse onto the angular momentum state 
— y + 2)- It is possible at this point to convince oneself that, iterating /c-times our procedure, 
we generate the coupled angular momentum state | y, — y + k). As far as the probability of 
success Pfe that after k measurements the N uncoupled spins are left in the state | ^ , — ^ + k) , 
it is easy to prove that it is given by 

k 

P k = ]Jsm 2 ( Pi at) (31) 

i=l 



where Vl = (f + 1) - (-f + i - l)(-f + i). 

Thus, if at the i— th step we have the possibility of choosing the time instant at which 

performing the measurement act on the central spin A in such a way that tn = ^2ap^ ' a ^ 

least in principle we may claim that the desired state | ^ , — ^ + k) of the N spins are generated 
with certainty. 



5 Estimating the order of magnitude of the coupling constant n ; 

As we are going to prove the dynamics of our system can be also successfully exploited in order 
to estimate the coupling intensity between each of the N distinguishable and non interacting 
N spins and the central one. 

Let's indeed consider for simplicity the case A = and, under this condition, concentrate 
on the behaviour of the probability Po(t) to find the system, at a generic time instant t, in the 
initial condition |-0o(O))- The results obtained in the previous Section immediately implies that 
such a probability is given by 



P (i) = |a(i)| 2 =cos 



2 



TV 



£°H- (32) 

i=i 



Thus, the probability of recovering the N uncoupled spins as well as the central one in their 
initial state is a periodic function of t, the period being in inverse relation to the quantity 
2j=i a | ■ 1^ as reasonable, we assume that all the coupling constants are of the same order, 
analyzing the temporal behaviour of Po(t), we have the possibility of estimating the order of 
magnitude of each otj . It is important to stress that the possibility of knowing at least the order 
of magnitude of the coupling constants aj play a central role for example in all the cases in 
which the spin of an electron localized in a quantum dot is used as realization of a quantum 
bit [22] . In these cases indeed the spin relaxation mechanism is mainly connected with its 
interaction with bulk nuclear spins. 

Let's moreover observe that the knowledge of the frequency of the function Po(t) given by eq. 
(|32p can be also exploited in order to estimate how much disordered the spin star system model 
(HJ) is. Let's suppose indeed that the N spins of interest have been prepared in a W-like state 
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following the procedure previously discussed. At this point, if we measure the observable a 3 z 
finding the eigenvalue +1, then the total system is projected onto the state |J., J., . . . , |j, • • • , I) 
in which only the j-th spin is in its upper state whereas the others are in their respective 
down state. The probability of such an event, directly obtainable starting from eq. (|14p . exactly 

a 2 

coincides with the quantity „ w 1 • 

> Q 2 

We may thus conclude that knowing the probability of success to find the j-th spin in the 
up state 1 1) when the system of N spins is prepared in a W-like state, allows us to give an 
estimation of the interaction strength between the j-th spin and the central one. 



6 Conclusion 

Generally speaking the possibility of establishing on demand fixed entanglement conditions 
in a multipartite system is an interesting objective both in its own and also in view of its 
applicative potentialities. In this paper in particular we have concentrated on a multipartite 
system composed by N not interacting spins |. In order to guide this system toward assigned 
entangled states, we have exploited the interaction between each of the N subsystems with 
a single spin i. The disordered spin star system thus obtained has been successfully used to 
generate W-like states as well as well-defined angular momentum states of the N uncoupled 
spins. 

The study of the exact dynamics of the disordered spin star system reported in this paper, 
has provided the possibility to envisage a way to estimate at least the order of magnitude of 
the coupling strength between each of the N uncoupled spins and the central one. To gain 
this information is, for example, of particular relevance when electron spin relaxation plays 
an important role. It is indeed appropriate to remark that our system can be adopted to 
describe hyperfine interaction of a single electron spin with nuclei in quantum dots and that 
this interaction mechanism may be the dominant source of electron spin relaxation. 
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